Fast ground-state cooling of mechanical resonator with time-dependent optical cavities 
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We propose a feasible scheme to cool down a mechanical resonator (MR) in a three-mirror cavity optome- 
chanical system with controllable external optical drives. Under the Born-Oppenheimer (BO) approximation, 
the whole dynamics of the mechanical resonator and cavities is reduced to that of a time-dependent harmonic 
oscillator, whose effective frequency can be controlled through the optical driving fields. The fast cooling of the 
MR can be realized by controlling the amplitude of the optical drives. Significantly, we further show that the 
ground-state cooling may be achieved via the three-mirror cavity optomechanical system without the resolved 
sideband condition. 

PACS numbers: 42.50. Vk, 03.65.-w, 42.50.Dv 
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I. INTRODUCTION 

Ground-state cooling of nano-mechanical resonators (MRs) 
has attracted great interests, as it is crucial in the improvement 
of detection precision of MRs llU2[], observation of quantum 
behaviors of macroscopic MRs [3], and quantum information 
processing based on MRs [4]. Over the years, methods to 
optimize the cooling of MRs in optomechanical systems [5- 
Fiql (or in electromechanical systems lfl7lo have been studied 
extensively. In these cooling schemes, the MR is coupled to 
a driven cavity through the optical radiation pressure and is 
cooled via the passive backaction cooling (also called self- 
cooling) JH> 0. ElrliHl without a feedback loop. More signif- 
icantly, the backaction ground-state cooling of MR lfl3l4l5ll . 
where the mean thermal occupation number of phonons n is 
less than 1, can be achieved when the resolved sideband con- 
dition is satisfied. 

Several experiments have reported progress in the backac- 
tion cooling of MR in optomechanical systems 151 4 1 ill , in par- 
ticular those approaching ground states of MRs. Recent re- 
search shows that MRs can be cooled to states with a mean 
thermal occupation number n = 35 in an optical cavity Holl . 
and n = 3.8 in a superconducting transmission line res- 
onator n liti . 

Most existing models for ground-state cooling of MR in 
optomechanical systems use resolved sideband cooling. The 
thermal phonon number of the MR is reduced by the inter- 
action between the MR and high-frequency auxiliary systems 
until the whole system eventually reaches equilibrium. Here 
we propose an alternative but efficient cooling method for MR 
in a three-mirror cavity (3MC) optomechanical system Il8ll . 
We notice that the MR and its auxiliary system, e.g., optical 
cavities, can be treated separately by the Born-Oppenheimer 
(BO) approximation [19]. The MR therefore behaves as a 
single-mode harmonic oscillator whose effective frequency is 
determined by the external optical driving fields. By control- 
ling the amplitude of optical driving fields, we can quickly re- 
duce the frequency of the bare MR, e.g., initially in a thermal 



equilibrium state, to a smaller effective frequency, but retain- 
ing the populations. In other words, the MR is cooled in a 
shorter time by doing work on the external system. Remark- 
ably, by combining such a fast cooling scheme with another 
process, we find that it is feasible to achieve the ground-state 
cooling of the MR with final effective frequency that is the 
same as the bare one. 

Our ground-state cooling mechanism is distinctly different 
from that of the conventional sideband cooling in optome- 
chanical systems: (i) Realization of sideband cooling requires 
a long time until a steady state is reached. Here the MR may 
be cooled down to a non-steady state at a short time and may 
later become hotter again because of interactions with envi- 
ronment, (ii) The optical driving field is time-independent in 
the sideband cooling, while the present scheme needs a time- 
dependent and much stronger optical power, (iii) The side- 
band cooling of MR happens when the optical detuning A ap- 
proaches the bare frequency of MR u, the cooling condition 
here is far detuned from the cavity (| A| 3> u>), which greatly 
simplifies the requirement for experimental control of the op- 
tical detuning, (iv) The decay rate of cavities in our scheme 
can be larger than the frequency of MR, while it is not allowed 
in resolved sideband cooling. 



II. THREE-MIRROR CAVITY CONFIGURATION 

Unlike the conventional two-mirror cavity (2MC) optome- 
chanical system, we consider here an optomechanical system 
of 3MC configuration, as shown in Fig.Q] The mirror of MR 
with two perfectly reflecting surfaces is placed inside a cavity 
and the two other fixed mirrors are transmissive and subject 
to external optical drives. This setup is different from that in 
Ref. H, where the MR is transmissive. 

The Hamiltonian of 3MC system can be written as 



H = H c + H s 



(1) 



where H s — p /2m + mwV/2 is the Hamiltonian for a MR 
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FIG. 1: (Color online) The schematics of 3MC optomechanical sys- 
tem. The movable MR with two perfectly reflecting surfaces is 
placed inside a driven cavity with two transmissive fixed mirrors. 



with effective mass m and bare eigenfrequency ui, and 

H c = ti(u) a — G a x)a^a + ti(u>b + Gbx)b^b 

+h{i a e l ^ t a + i b e tVbt b + H.c.) (2) 

describes the optomechanical coupling between single-mode 
optical driven cavities a and b (with the eigenfrequencies u> a 
and Ldf,, respectively) and the MR. Here £ a j, are the ampli- 
tudes of external optical driving fields to the cavities with 
the corresponding drive frequencies v a ^, G a = u) a jL a and 
Gb = ^>b /Lb are the corresponding optomechanical coupling 
strengths via the radiation pressure, and L a .b are the lengths 
of the cavities a and b. 

Similar to the case for electrons and nuclei in a molecule, 
the eigenfrequencies u! a and LJb of two cavities are much 
higher than the frequency lu of the MR such that the BO ap- 
proximation can be employed to separate the degrees of free- 
dom of MR from those of two cavities. For simplicity, but 
without loss of generality, we set u a = Ub, G a = Gb = G, 
L a = L b = L, £ a = £* = £ b = 4* = 6 and v a = v h = v. 
The cavity Hamiltonian can now be written as 

H c = h{A-Gx)a''a + h(A + Gx)b''b + h((a + b+H.c.) (3) 

in a rotating frame with respect to Hf ,o — hv(a^a + b%), 
where A = (oj a — v) is the optical detuning at the absence of 
the motion of the MR. 

We define operators a = a — a and b = b — f3, where 
a = -£/(A - Gx) and /3 = -£/(A + Gx), such that the 
Hamiltonian @ can be expressed as 



H c = h(A - Gxpa + h(A + Gx)tfb + c (x), (4) 



where 



c (x) 



\ 2 - G 2 x 2 



(1 



A 2 L 2 ' 



is x dependent. Note that we have used the condition 
Gx/A <C 1 and ignored the effects from the decay rates 
(k„ = Kb = n) of cavity fields based on the fact that k is 
much less than the detuning A. 

When the cavity fields are fixed for the state |0) (= \0(x))) 
with a|0) =0 and 6|0) = 0, the effective Hamiltonian of the 
MR is 

(P + Arf 1 



H cS = " ~~ u; +-muj 2 x 2 - 
2m 2 

P 2 1 2 2 
2m 2 



c (x) 



where the effective frequency is given by u> 2 s — lj 2 — 
Am^u?J(mA z L 2 ) and c r = -2h 2 ^ 2 /A is a constant shift. 
We have assumed that the induced gauge potential = 

-^(0|V|0) = -ih(0(x)\J^\Q(x)) = if |0(.t)) is a real 
function. Then the MR behaves as a harmonic oscillator with 
the effective frequency determined by £ (equivalently the am- 
plitudes of the optical driving fields), which is similar to the 
discussion according to the conventional optical springs B20I1 . 

The 3MC configuration distinguishes itself from the 2MC 
configuration because the high symmetry of the three-mirror 
cavity allows to significantly improve the trap stiffness and 
to partially remove the bistability as addressed in Ref. lfl8ll . 
The MR displacement in the effective Hamiltonian (O is zero, 
while the corresponding MR displacement is finite in the 2MC 
configuration, making the dynamics of the 2MC more compli- 
cated. 



III. GENERAL FAST COOLING OF TIME-DEPENDENT 

MR 

Although the effective Hamiltonian (0 is obtained by us- 
ing the time-independent BO approximation for the time- 
independent £, the similar result also holds for a time- 
dependent parameter = £o/(i) when the factor f(t) 
(\f{t)\ < 1) varies very slowly in comparison with the de- 
tuning A, i.e., f{t)/f{t) -C A. In this case, the Schrodinger 
equation H l^} = ifidt \ip) can be simplified by a simi- 
lar BO approximation with a highly non-trivial considera- 
tion of time dependency, as discussed in Ref. [21]. The 
BO approximation implies (0\H s \ip) rs H s (0\ip) such that 
(0\H\ip) « (H s +cq(x))(0\i[i) . Here (0\ip) is an instantaneous 
eigen wave-function of the MR with the effective Hamiltonian 
H c , cq(x) is the corresponding instantaneous eigenvalues and 
now is time dependent, cq(.t) = cq(x, t). 

Based on the adiabatic theorem with conditions 
f(t)/f(t) <g; A and the wave function |0) being real, 
we obtain <9t(0| s» 0. The effective Schrodinger equation for 
the MR is then given by 



[H B +c (x,t)]@\il>) = ihd t @\il>). 



(6) 



Therefore, the MR can be described as a time-dependent har- 
monic oscillator with the effective Hamiltonian 



H eB (t) =H S + Co(x,t) = 7^ + ^muj 2 s (t)x 2 . 
Here the corresponding time-dependent eigen-frequency 
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Weff(t) = LO^/l + l]f 2 (t) 



(7) 



(8) 



(5) 



with r\ = —4h^ 2 u} 2 /(mijj 2 A 3 L 2 ), is controlled by the exter- 
nal optical field via the dimensionless function f(t). 

The dynamics of a time-dependent harmonic oscillator is 
analytically solvable using Lewis-Reisenfeld invariants l22tl . 
Here we are only interested in the specific trajectories u> e s(t) 
wherein the instantaneous populations at the initial time are 
the same as those at the final time, but the eigenfrequency and 
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the corresponding average energy at the final time decreases. 
That means the harmonic oscillator has the same entropy at 
the initial and final times and is cooled by doing work on the 
external field, instead of having heat flowing out of it. Such 
a cooling trajectory does not depend on the initial state. The 
simplest way to realize the cooling trajectory is the adiabatic 
process where the populations remain the same all the time. 
The disadvantage is that the adiabatic process needs a long 
time, during which the relaxation of the harmonic oscillator 
may bring negative effects. Alternatively, an optimal bang- 
bang process has been proposed to achieve a cooling trajec- 
tory in a shorter time [23]. The effective frequency is real and 
stepwise function of time in this process, where the final in- 
stantaneous populations of the MR remain the same as the ini- 
tial ones in terms of controlling the step values and durations 
of the effective frequency. Recently a fast optimal friction- 
less process l24ll25ll was proposed to achieve atom cooling by 
choosing a certain trajectory of effective frequency, wherein 
the instantaneous effective frequency is allowed to be "imagi- 
nary". 

Employing these methods, one is able to cool the MR in 
time-controllable optomechanical systems, as done in atom 
cooling 1I23TI25I1 . Let us estimate the feasibility of cooling a 
MR with reasonable empirical parameters: the bare eigenfre- 
quency of MR uj/2tt — 134 kHz, the effective mass m = 50 
Pg iHD, the eigenfrequency of the optical cavity u a /2Tr pa 
7 x 10 14 Hz, the optical detuning A/27T pa 10 7 Hz, and the 
length of the optical cavity L pa 2 mm. The optical driving 
fields are controllable and can have the time-dependent form 
£ = &/(*) with £ /27T w 10 9 Hz and < 1. With 

these parameters, we obtain 77 pa —9. The effective frequency 
w c ff (t) can be much smaller than the bare frequency u) of the 
MR, or even be "imaginary" by tuning /(£). 

Specifically, consider a MR in the thermal state 
with the mean phonon occupation number n(ij) = 
l/[exp(^w/fc B r) - 1] pa k B T/huj (> 1) at the initial time 

and cj e s(ti) = u. One can design a trajectory between 
u) e s(U) an d the much lower effective frequency u) e g(t{) at 
the final time (where uj e g(t{)= oj/R with R 3> 1 and the 
corresponding /(if) = a/(1 — R~ 2 ) /rf) such that the final 
mean occupation number n(tf) = n{ti). The energy of the 
MR is decreased by a factor R after doing work on the ex- 
ternal field, and the final effective temperature is reduced to 

Teff(tf) = Mtf)Weff(tf)/fcB = T I R « T - 

We have proposed a fast cooling mechanism for the MR as 
a time-dependent effective harmonic oscillator in a 3MC op- 
tomechanical system, where the final effective frequency of 
the MR is different from the initial one. However, the ap- 
proach will not lead to the ground-state cooling since the final 
mean population number of the effective harmonic oscillator 
is the same as the original one which is usually much larger 
than 1. 



IV. GROUND-STATE COOLING OF MR 

In the preceding section, we obtain a fast cooling of MR 
wherein the final effective frequency of the MR is different 



from the initial one. People may be mostly interested in the 
cooling of the MR with the final effective frequency of the 
MR being (nearly) equal to the original bare one, such as the 
sideband cooling of MR. In this section, we design a different 
cooling approach for the same 3MC optomechanical system to 
cool the MR close to its ground state in a short time, with the 
final effective frequency of the MR being equal to the original 
bare one. 

We consider that the MR is originally (at 
t = t Q ) in a thermal equilibrium state p(t ) — 
e -H c!! (t )/k B T / Tr ( e -ffcf f (to)/fc B T^ i n a bath at a refrigeration 

temperature T = 20 mK [26] in the absence of the optical 
driving fields. The corresponding mean occupation number 
of thermal phonons is n(t Q ) = l/[cxp(huj/kBT) — 1] ^> 1. 
The MR is then subject to the optical radiation pressure of 
the driven cavities, described by the effective Hamiltonian 
in Eq. Q. We use the same preceding typical param- 
eters, except with the optical detuning A/2n = — 10 7 
Hz and with the maximum value of optical driving fields 
£(°)/2tt pa 10 12 Hz, and obtain 77 « 9 x 10 6 . At time 

t = t\ the MR evolves to a new thermal state p(t\) — 
e -H ea (ti)/k B T ^ TT ^ e -ff etf (ti)/fc s T') with the mean dermal 

occupation number n(ti) = l/[exp(Uu} e g(ti)/ ksT) — 1]. 
We take f(t{) = 1 such that u; e ff(*i) ~ 3000w and 
n(t{) pa 0.66 < 1 <C n(t Q ) «a 3200. The intermediate 
trajectory between t a and ti may be arbitrary, on the condition 
that the MR is in the thermal state at time ti. 

We need to design a special trajectory of the effective fre- 
quency to e ff(t) from t = ti to the final time t = tf such that 
w e g (tf ) equals to the original bare frequency oj and the corre- 
sponding state p(tf) — p(ti). Therefore, the final mean ther- 
mal occupation number n(tf) = n(ti) pa 0.66, which indi- 
cates that the MR is cooled near its ground state. 

We now focus on designing such a trajectory of w e g (t) 
in terms of the Lewis-Riesenfeld invariant of motion 12211 . 
The invariant of motion for our MR harmonic oscillator is 
I(t) = mu} 2 x 2 /[2b 2 (t)} + [b(t)p - mb(t)x} 2 /2, where u = 
ui e s(U). The dimensional real function b(t) satisfies the con- 
dition b(t) + u 2 s (t)b = oj 2 /b 3 (t) JH, which is used in 
Ref. 1I24I1 for cold atoms. We show that the simplest polyno- 
mial choice of b(t) 112711 will lead to an above-required cooling 
trajectory, which can be realized by controlling the parameters 
f(t), as shown by the solid and dashed lines in Fig. [2] 

In addition, the cooling process could even be faster if 
imaginary frequencies w c g(t) could be allowed as done in 
Ref. 0241], where the imaginary frequencies cooling scheme 
were first proposed. The single optical mode Hamiltonian (O 
fails to provide such imaginary values of u! e g(t) since 77 is 
positive for the fixed parameters above. For that purpose, we 
could adopt two optical modes in each cavity for the simi- 
lar symmetric 3MC configuration. A self-adjoint Hamilto- 
nian B28I1 for the model with two optical modes in each cavity 
in a rotating frame should be 

H' c = h{At -Gxx)a\a x + ft(Ai +Gix)6 t 1 fei 
+?i(A 2 - G 2 x)ala 2 + h(A 2 + G 2 x)b\b 2 
+fc[£i(ai +61) +6(02 + 6 2 )H.c], (9) 
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where ai.2 (^1,2) denote the annihilation operators for the two 
modes in cavity a (cavity b) with the frequencies u)\ <z. The 
coupling strengths of optical radiation pressure are G12 = 

u>i,2/L. :=^ 0) /i(t) and 6 := ^ f 2 (t) (|/ 1)2 (t)| < 1) 
are the amplitudes of the external time-dependent optical driv- 
ing fields with frequencies v\fl. The optical detunings are 

Ai = (uji - vi) and A 2 = [u% - v 2 ). 

Following the preceding discussion, the effective time- 
dependent Hamiltonian for the MR becomes 



P 

2m 



1 



o? S (t)x 2 , 



where the time-dependent eigenfrequency is given by 



u>3?(t) = " 2 li + m fi(t) + mfi(t)] 



(10) 



(ii) 



with two coefficients 771 = —4H^[ lo 2 / {muj 2 /S.\L 2 ) and 
m = -4^ 0)2 w|/(mw 2 A3i 2 ). 

We now have more effective control over the time- 
dependent eigenfrequency w' eS (t) through two independent 
parameters /i,2(t) (proportional to the external input optical 
powers) and will use the same typical parameters for both 
modes as before, except for the optical frequencies lo\/2-k « 
7 x 10 14 Hz, 0J2/2-K » 6 x 10 14 Hz, the maximum amplitude 



of optical driving fields ^ 0) /2tt w 10 12 Hz, /2ir w 10 10 
Hz, and the detunings Ai/27r = -10 7 Hz and A 2 /27r = 10 7 
Hz. Therefore, the coefficients iji «9x 10 6 and r/2 « —670. 
Thus, the effective frequency could be imaginary in the pro- 
cess when f% (i) is much larger than f\ (t) and the cooling can 
be achieved even faster (see the blue dotted and dashed-dotted 
lines in Fig.[2)- 

We wish to mention that our ground-state cooling scheme 
requires a strong input optical power to obtain a strong modu- 
lation for the effective frequency. The maximum value of the 
optical input power is P = ?kj a £(°' 2 /2k ~ 1 W for the max- 
imum £(°) /2tt w 10 12 Hz, where the decay rate of the optical 
modes of the cavities is assumed as n/2n w 10 6 Hz. The 
corresponding high finesse cavities with F = ttc/(2Lk) « 
3.8 x 10 4 , are required. We remark that the adiabatic condition 
for the time-dependent BO approximation , fit)/ fit) <C | A| 
an d fxpit)/ fiflif) *C IA12I, is satisfied for the trajectories 
in Fig. |2] 

Our general analysis shows that we can realize the ground- 
state cooling for a MR with the present approach. Even if the 
coupling between the MR and the optical cavity is removed 
by shutting off the optical driving fields, the MR will remain 
in the same ground state for a longer time after t{, until it 
approaches to a new thermal equilibrium with its bath. Dur- 
ing this period of time, one can perform quantum operations, 
such as quantum information processing, on the MR. Note 
that the controlled trajectory between uj s s(ti) and ui s(t{) is 
preformed in a short time of a fraction of 1/lu, such that the 
relaxation process of the MR with a high Q factor (e.g., 10 5 ) 
can be ignored. 
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FIG. 2: (Color online) (a) The square of instantaneous eigen- 
frequency w' B ff(t) (or w e ff(£)) of MR from u/ eff (£j = 0) = 3000oj 
to Weff(if) = u) for different durations: tf — 0.2 /lj (blue dotted 
line); tf — 0.6 /lj (black solid line); tf = 2/lj (red dashed line). 
Here, the bare frequency of MR is lj = 2ir x 134 kHz. The nega- 
tive value is due to an imaginary w eff (t) which appears for the case 
of short evolution time tf. (b) The corresponding time-dependent 
control parameters / 2 (t) [flit)] (proportional to the optical powers) 
and /| it) (the blue dotted-dashed line, which contributes mainly to 
the "imaginary" u}' eif it)). 



V. CONCLUSIONS 

We propose a feasible scheme to cool a MR through a 3MC 
optomechanical system. In the limit of large optical detuning, 
we show that the optical fields can be eliminated and the de- 
gree of freedom of the MR in the general Hamiltonian can 
be described by a time-dependent harmonic oscillator. By 
controlling the amplitude of the external input optical driv- 
ing fields, we can manipulate the time-dependent effective 
frequency and obtain the cooling of the MR via keeping the 
same populations of instantaneous levels at the initial and fi- 
nal times during a special trajectory of short time period. The 
ground-state cooling of MR can be realized as well. It is en- 
couraging that a similar trajectory has been implemented ex- 
perimentally M25I1 for atom cooling. We believe that the same 
technique is applicable to our scheme for a mechanical res- 
onator. 



Acknowledgments 

We would like to thank Dr. P, Zhang for helpful discussions. 
This work was supported by the RGC grant of Hong Kong 
No. HKU7044/08P, the Ikerbasque Foundation Start-up, the 



5 



Spanish MEC No. FIS2009-12773-C02-02, and the State Key Program for Basic Research of China (No. 2006CB921800). 



[1] V. B. Braginsky and A. B. Manukin, Measurement of Weak 
Forces in Physics Experiments (The University of Chicago 
Press, Chicago, 1977). 

[2] C. M. Caves et al., Rev. Mod. Phys. 52, 341 (1980); M. F. 
Bocko etal, ibid. 68, 755 (1996); M. D. LaHaye etal, Science 
304, 74 (2004); B. Abbott et al, Phys. Rev. Lett. 95, 221101 
(2005). 

[3] S. Mancini et al, Phys. Rev. Lett. 88, 120401 (2002); W. Mar- 
shall et al, ibid. 91, 130401 (2003); L. F. Wei et al, ibid. 97, 
237201 (2006); F. Xue et al, Phys. Rev. B 75, 033407 (2007). 

[4] S. mancini et al, Phys. Rev. Lett. 90, 137901 (2003); K. Ham- 
merer et al, ibid. 102, 020501 (2009); L. Tian and P. Zoller, 
ibid. 93, 266403 (2004). 

[5] C. H. Metzger and K. Karrai, Nature 432, 1002 (2004). 

[6] S. Gigan et al, Nature 444, 67 (2006). 

[7] O. Arcizet et al, Nature 444, 71 (2006). 

[8] J. D. Thompson et al, Nature 452, 72 (2008); Sankey et al. 

Nature Physics 6, 707 (2010). 
[9] A. Schliesser et al, Nature Physics 4, 415 (2008). 
[10] S. Groblacher et al, Nature Physics 5, 485 (2009). 
[11] T. Rocheleau et al, Nature 463, 72 (2010). 
[12] S. Mancini et al, Phys. Rev. Lett. 80, 688 (1998). 
[13] I. Wilson-Rae, N. Nooshi, W. Zwerger, and T. J. Kippenberg, 

Phys. Rev. Lett. 99, 093901 (2007). 
[14] F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, Phys. 

Rev. Lett. 99, 093902 (2007). 
[15] Y. Li, Y. D. Wang, F. Xue, and C. Bruder, Phys. Rev. B 78, 

134301 (2008); Y. Li, Z. D. Wang, and C. P. Sun, Eur. Phys. J. 

D61, 215 (2011). 
[16] T. J. Kippenberg etal, Optics Express 15, 17172 (2007). 
[17] L. Tian et al, Phys. Rev. Lett. 93, 266403 (2004); P. Zhang 

et al, ibid. 95, 097204 (2005); A. Naik et al, Nature 443, 

193 (2006); Y.-D. Wang et al, Phys. Rev. B 80, 144508 

(2009); K. Jacobs, H. I. Nurdin, F. W. Strauch, and M. James, 

larXiv: 1003.26531 



[18] M. Bhattacharya and P. Meystre, Phys. Rev. Lett. 99, 073601 
(2007); P. Meystre et al, J. Opt. Soc. Am. B 2, 1830 (1985). 

[19] M. Bom and J. R. Oppenheimer, Ann. Physik 84, 457 (1927). 

[20] B. S. Sheard et al. Phys. Rev. A 69, 051801(R) (2004), and 
references therein. 

[21] T. J. Gil et al, Phys. Rev. A 45, 6479 (1992). 

[22] H. R. Lewis et al, J. Math. Phys. (N.Y.) 10, 1458 (1969); M. V. 
Berry et al, J. Phys. A 17, 1805 (1984); V. V. Dodonov et al, 
J. Math. Phys. (N.Y.) 34, 3391 (1993); M. A. Lohe, J. Phys. A 
42, 035307 (2009). 

[23] P. Salamon etal, Phys. Chem. Chem. Phys. 11, 1027 (2009). 

[24] X. Chen et al, Phys. Rev. Lett. 104, 063002 (2010); J. G. Muga 
etal.,]. Phys. B 43, 1 (2010). 

[25] J.-F. Schaff etal, Phys. Rev. A 82, 033430 (2010). 

[26] We notice that the conditions of initial temperature 20 mk and 
the strong optical input power 1 W may be at or beyond the 
borderline of current optomechanical experiments. We choose 
such ideal parameters to show how the ground-state cooling of 
the MR with bare frequency u — 2tt x 134 kHz can be reached. 
For the MRs with larger bare frequency and smaller effective 
masses, the less laser power and higher initial temperature are 
required to achieve its ground-state cooling. 

[27] That is, b(t) = 6C[t/( tf-ti)] 5 -15C[£/ (>-ti)] 4 + 10<[t/(tf- 
ti)] 3 + 1 with ( = V^fi(ii)M=ff(tf) - !■ 

[28] We emphasize here that the Hamiltonian is self-adjoint when 
the "eigenfrequency" of a harmonic oscillator is imaginary. 
While we do this as an additional exercise for comparison, our 
MR cooling scheme works well without the requirement for 
imaginary eigenfrequencies, as shown in Fig. [2] (the black solid 
and red dashed lines), which will be for the practical implemen- 
tation of ground state cooling of MRs and whose rate is still 
comparable with the cooling based on the bang-bang process as 
in Ref. (H. 



